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Abstract 

An a-helical region of protein molecule was considered in a model of nanotube. The molecule is in conditions of 
quantum excitations. Such model corresponds to a one-dimensional molecular nanocrystal with three molecules in 
an elementary cell at the presence of excitation. For the analysis of different types of conformational response of 
the a-helical area of the protein molecule on excitation, the nonlinear response of this area to the intramolecular 
quantum excitation caused by hydrolysis of adenosine triphosphate (ATP) is taken into account. It has been established 
that in the simplest case, three types of excitation are realized. As estimates show, each of them 'serves' different kinds 
of protein. The symmetrical type of excitation, most likely, is realized in the reduction of traversal-striped skeletal 
muscles. It has the highest excitation energy. This well protects from casual actions. Antisymmetric excitations 
have intermediate energy (between symmetrical and asymmetrical). They, most likely, are realized in membranous 
and nucleic proteins. It is shown that the conformational response of the a-helical region of the protein is (in 
angstroms) a quantity of order N c /5, where N c is the number of spiral turns. For the number of turns typical in this 
case: A/ c ~ 10, displacement compounds are a quantity of order 2 A. It qualitatively corresponds to observable 
values. Asymmetrical excitations have the lowest energy. Therefore, most likely, they are realized in enzymatic 
proteins. It was shown that at this type of excitation, the bending of the a-helix is formally directed to the opposite side 
with respect to the antisymmetric excitations. Also, it has a greater value than the antisymmetric case for N c < 14 and 
smaller for N c > 14. 
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Background 

Hydrolysis of ATP and amide I excitation 

A protein molecule has a rather unique structure not only 
in the chemical-biological point of view but also as an in- 
teresting physical and mathematical object. If we consider 
it as a physical object, then such object may be referred to 
as a nanostructure without any doubt. Thus, the alpha- 
helical region of a protein molecule simultaneously may 
be considered both as a nanotube and as a nanowire: this 
depends on the considered level of structure. 

Here, the alpha-helix is considered at the level of second- 
ary structure where it is a nanotube. It is in the conditions 
of quantum excitation which is stimulated by reaction of 
hydrolysis of adenosine triphosphate (ATP). As a result of 
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this reaction, energy in the form of quanta of infrared range 
is released. It is considered that they are absorbed by a 
group of energy states known in an alpha-helix as amide I, 
etc. It is considered also that these absorbing states have an 
internally molecular oscillating nature. The results obtained 
here allow giving a definite answer to this question, because 
in the infrared range, absorption can also have the nature 
of electronic transitions between states with the main 
quantum number equal to 2. 

The alpha-helix is interesting as a mathematical object 
too. Due to the high sensitivity of its crystalline lattice' 
in relation to excitation, we are coming to a necessity to 
solve a nonlinear system of the so-called eigen type, i.e., 
actually, we are coming to a necessity to search for the 
eigenvalues and eigenvectors of a nonlinear system of al- 
gebraic equations. Such a problem, as it is known to us, 
is a scantily explored mathematical problem. 



Spri 
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Figure 1 shows the alpha-helical fragment of a protein mol- 
ecule. Similar regions in proteins are widespread enough 
in vivo. The degree of helicity in different proteins varies from 
12% to 96%. As can be seen from Figure 1, the alpha-helical 
fragment of protein molecules is structurally a nanotube. The 
same is true for its physical properties. Therefore, to such re- 
gions of protein molecules in their excited states, it is natural 
to apply methods that are specific for nanotubes. 

As a result of hydrolysis of ATP molecule, energy is 
realized in the range 0.2 to 0.4 eV\ It depends on the 
charge state of the ATP molecule, in which the compos- 
ition of the environment influences mainly (pH, etc.). 
The energy of hydrolysis is absorbed by an alpha-helical 
region of the protein molecule. It takes place due to in- 
ternal vibrational excitations of the peptide groups 
(HCNO) in the state amide I. Its energy is also varied 
within the limits of 0.2 to 0.4 eV. These excitations in- 
duce a significant increase of dipole moments of the 
peptide groups, which is equal to 3.7 D, on 0.29 D [4,5]. 

There exists another point of view. Excitation of amide I 
may have an electronic nature. It may correspond to tran- 
sitions between energy bands with principal quantum 




Figure 1 The real (a) [1,2] and schematic (b) [3] images of an 
alpha-helix. 
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numbers that are equal to 2. The physical nature of excita- 
tion is inessential for further calculations, but further it 
will be shown that their nature may be determined 
experimentally. 

Methods 

Amide I excitation in the simplest model of alpha-helical 
region of protein 

Foremost, we need to determine the model of descrip- 
tion of the spatial structure of the alpha-helix. Since it is 
considered as a molecular crystal, the nearest neighbor 
approximation is used, which is typical for such crystals. 
However, as seen from Figure lb, the nearest neighbors 
for some peptide group with number n are not only 
group n ± 1 but also group n ± 3. 

The simplest model of the spatial structure of the 
alpha-helix is shown in Figure 2. Such simplified model 
differs from a real molecule only by symmetry. In the 
model considered, the molecule is independent from 
each other: translational and axial symmetries. The real 
molecule has translational-helical symmetry. Preliminary 
investigations have already shown that the qualitative 
picture in terms of types of excitation does not change. 
Changes will only be quantitative. They will lead only to 
some displacement of the absorbing states by energy 
and, may be, to some mixing of states. In this, simplest, 
model, all turns of the helix closed on itself, although 
Figure 1 shows that this is not quite so. Each turn of the 
helix is open for the nearest neighbor. It was previously 
shown [6] that taking into account open individual cells 
leads only to quantitative changes. The qualitative pic- 
ture remains unchanged. 

Within the framework of the considered model, every 
three peptide groups that belong to one turn of the helix 
grouped into one complex unit cell. We will number 
these unit cells by indices n, m, etc. The number of such 
cells is three times less than the number of peptide 
groups, i.e., N 0 /3. Peptide groups within a single cell will 
be enumerated by indices a, ft, etc. that may take values 
0, 1, 2. The general functional for the alpha-helix in this 
model has the form [7] 




Figure 2 Simplest model of alpha-helix as a one-dimensional 
molecular crystal with three molecules per unit cell. Arrows are 
showing a separate peptide group. They symbolize the 
dipole moments. 
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w(R na - R m ^) in this functional is the basic energy of 
interaction between peptide groups na and m/?. It is in- 
dependent on the presence of excitation and exists al- 
ways. D(R na -R m p)\A an \ 2 is an additional energy to the 
w(R na - R m p) energy of interaction related only to excita- 
tion but considerably smaller. Factor A an is the wave 
function that describes the excited state of the examined 
alpha-helical region of the protein molecule. It deter- 
mines the spatial-temporal distribution of excitation in 
this region. The energy D(R na - R m p)\A an \ 2 leads to the 
breaking of the equilibrium of the alpha-helix and stimu- 
lates its conformational response to excitement. Energy 
Mna^pApnjAan is also an additional energy of interaction. 
However, it is much less than D(R na - R m p)\A an \ 2 but 
important because it provides the propagation and trans- 
fer of excitation along the alpha-helix. 

As shown in Figure 2, the nearest neighbors for some 
peptide group na will only be the peptide groups m = n±l, 
fi = a and rn = n, f> = a±l. Taking into account that in the 
considered model all energy terms depend on the distances 
between amino acid residues only, the following formulae 
in the nearest neighbor approximation may be obtained: 



v w + l,a 



R, 



n,a\) rna 



Pna — | R«,a + 1 ^n,a I ♦ 



Rnc 

Let us take into account that the response of the lat- 
tice (Figure 2) on excitation inside of the unit cell is 
small enough. Thus, it may be neglected in comparison 
with a similar response between unit cells. In this sense, 
the equality p na = p 0 is always supposed fulfilled. Factor 
R na is the only value that takes into account the response 
of the alpha-helix on excitation. Thus, we will denote its 
equilibrium value as R 0 . Values p 0 and R 0 are shown in 
Figure 2. Taking into account the normalization condition 



Ew 2 = i, 

na 

the last functional takes the form 



(i) 



E({R},{A}) = +£>± + { W ( R na)+D(Rna)\A a „\ 2 + 

na na 

+ +<,n-l) A an+-M ± (A* a+ln +A* a _ 1 ^A an ) • 

(2) 

Here, w ± = w(p 0 ), D ± = D(p 0 ), M ± = M(p 0 ), and M\\ =M 
(R 0 ). Obviously, \M ± \ *|M|||. If resonance interaction has 
no electronic nature, inequality will be realized: \M ± \ < \M\\\. 
If excitation has an electronic nature, inequality will be 
reversed: \M ± \ > \M\\\. This difference may be detected 
experimentally, and the answer of the question about 
the physical nature of excitation may be obtained. 



New equilibrium values of distances, which actually 
coincide with the step of alpha-helices, are determined 
using the general condition of minimization: dE ^^ A ^ = o. 
When interactions between peptide groups are modeled 
as purely dipole, the step of the alpha-helix always de- 
creases and is given by 



\D/{Ko)\ 

Kna — , /// n X, \^na\ • 

\w//(Rq)\ 



(3) 



Next, we must substitute (3) in (2), take into account 
the condition j^/^J \A na \ 2 « R 0 , designate w(R 0 ) = w\\, 



=G , and introduce convenient re- 



designation: M\\ = -\M\\\ = -2A, M ± = \M ± \ = 277, which 
take into account the true signs. Then for the functional 
(2), finally, the following formula will be obtained: 

E({A}) =E om -J2 {a(K,„+x +K,n-iy m +\G\A an \ i 

na 

- n ( A : + l,n+ A U,n) A an}- 

(4) 



In Equation 4, E OCii = (w ± + W\\)N 0 + D ± + D\\, and the 
following is taken into account: 

+ = (w ± + W\\)Nq. 

na 

N 0 is the number of amino acid residues in the alpha- 
helical region of the protein molecule, which is under 
consideration. 

Further, for implementation of the conditional minimization 
of energy (4) in relation to wave functions A aw it is ne- 
cessary to create a conditional functional: E yM ({A}) = E 

({A}) + s\ |A^| 2 ♦ From a mathematical point 

\ na J 

of view, parameter s is an indefinite Lagrange multi- 
plier, and physically, it is the eigenvalue of the consid- 
ered system. The minimization procedure dE * d }* A ^ = 0 

produces the equation A(A a>n + 1 + A a>n _ i) + G\A an \ 2 
A an - n(A a + hn +A a _ hn ) + sA an = 0. After dividing this 
equation by A and introducing the notations, 



77 



A 1 



(5) 



it is possible to reduce it to a dimensionless form: 



A a ,n+1 ~\~A a ^ n -\ + g\A an | A an -XA a^i^-XA a -l,n ^X^-an — 0. 

(6) 
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The function A an is complex. Therefore, the common 
solution of the system (6) has the form A an = a an • exp 
(iy an ). Amplitude a an and phase y an are real functions of 
the variables a and n. We confine ourselves to the 
Hamiltonian-Lagrangian approximation in phase [8]. 
Due to the stationarity of the solved problem, this ap- 
proximation has the simplest form: y an = kn. If the 
alpha-helical part of the molecule is long enough, b a 
Born-Karman condition gives k — ^j. Here, N c =^- is 
the number of turns in the considered alpha-helical re- 
gion of the protein molecule. It plays the role of the di- 
mensionless length of the helical region of the protein 
in units of an alpha-helix step. Parameter / has the 
values ;«0, ± 1, ± 2, Then 



Jkn 



(7) 



and Equation 6 takes the form 



-ie 



Separating real and imaginary parts, we have the fol- 
lowing formulae: 



cos (k)'( 



sm(k)-(a a , n+ i-a a ,„-i) = 0. 



(8) 
(9) 



The solution of this system is usually determined after 
transition to continuous approximation. But we will 
analyze systems (8) and (9) without using the continuous 
approximation, because we are interested in very short 
alpha-helical regions (10 to 30 turns). 

There is only condition a a>n + ± - a a , n _ 1 = 0 (if not to 
restrict solutions by using the condition /c = 0), which 
does not depend on any symmetry of the alpha-helix: 
whether it is the symmetry of the model or the sym- 
metry of the real molecule. Viewing of other conditions 
can appear useful on account of the real structure of the 
alpha-helical region. In the simplest case, it may be re- 
duced to the equation a an = P a . The system (8) now de- 
generates in the system of three nonlinear equations: 



x p 0 -P 1 -P 2 +yH = 
-P 0 -P 1 +xP 2 



0; 



-yPl 



0; 
0; 



(10) 



o ' 



1 



where the following designations are introduced: 



j + 2cos(/c) 



g. 
A" 



(11) 



The last, fourth, equation arose out from normalization 
condition (1). The coefficients P a {a = 0, 1, 2) determine 
the excitement of each peptide chain as a whole. 

The system (10) consists of four nonlinear equations 
for determining the values P 0 , P 1} and P 2 and the eigen- 
value x. By adding and subtracting the first two equa- 
tions and some transformation of the third equation, the 
system (10) can be reduced to the form 



(p 0 -P 1 )[x+l+y(P 2 0 +P 2 1 ArP 0 Pi)] =0; 

(p 0 + y(Pl +Pj-PoPi)} = 2P 2 ; 

(P 0 +i>i) = (x + yPi)P 2 \ 



(12) 



Pi 



Pl+P 2 2 



1 



This transformation does not affect the solutions of 
the system. 

For the solution, the condition P 0 + P 1 = 0 should be 
used. This condition together with the condition P 2 = 0 
turns into an identity the second and third equations. 
After some simple transformations, we obtain the anti- 
symmetric excitations: 

n(«) _ 1 . nW _ 1 . 



Ja) 



0; 



-1- 



2N C 



Using Equations 4, 5, and 11, it is possible to find the 
energy: 

G 



E a (k) 



+ e a {k) =E QI 



-77- 



2N C 



-2Acos(k). (13) 



Next, we use the condition P 0 -P 1 = 0, which turns 
into an identity the first equation in (12). After some 
analysis, we can find two types of excitation: 

• Symmetrical 



p( c ) _ p( c ) _ p( c ) 

1 0 — 1 l — 1 2 



3AL 



For these excitations, in analogy to the 
antisymmetric, it is possible to obtain the energy: 



E c (k) = £ 0CH + 2/7- — -2A cos(/c). 
• Asymmetrical 



(14) 



2y_ 
3N C 



For these excitations, it is also possible to get energy: 
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E H (k) =E m -n- — -2Acos(k). (15) 



The energies E a (k), E c {k), and E H (k) contain parameters 
A = |Af|| |/2 and TI=\M ± \I2. As it was noted between 
Equations 2 and 3, the relation between these parame- 
ters makes the determination of the physical nature of 
excitation possible: whether they are electronic or intra- 
molecular. Because one of them (A) determines the 
width of the excited energy bands, and the other (77) 
their positions, this is the basis for the experimental ana- 
lysis of the nature of excitations. 

There are a few possibilities else for searching for solu- 
tions of the system (12). Preliminary analysis shows that 
the obtained excitations are peculiar in a more or less 
degree for both symmetries: whether it is the symmetry 
of the model or the symmetry of the real molecule. The 
other solutions of the system (12) need to be analyzed 
only in the conditions of the maximum account of the 
real structure of an alpha-helix. But the general analysis 
of this system shows that the solutions of a new quality 
are not present: all of them belong to the asymmetrical 
type. However, attention should be paid to the equation 
a a,n + i ~ a a,n -i = 0, which has led to the requirement 
Man = Pa- This condition is strong enough and essentially 
limits the solution: it is a constant in variable n, i.e., does 
not have the spatial distribution along an alpha-helix. 

Results and discussion 

The analysis of the energetics of the protein excitation 

From definitions (13), (14), and (15), it ensues that received 
excitations are located in accordance with the inequality E c 
(k)>E a (k)>E H (k). Thus, E c (k)-E a {k) = 3/7 + ^f- and E a 
(k)-E H (k) = ^r. It can be seen that for the alpha-helical 
region of finite length, when the number of turns N c * °o, 
the lowest energy is the energy of asymmetric excitation 
E H . Also, it is visible that energy E c is always strongly sepa- 
rated from energies E a and E H . Even when the number of 
turns N c ^oo and the energies E a and E H practically coin- 
cide, the energy E c is separated from E a and E H on a value 
3/7=3|Mj/2. Amide I excitations manifested experimen- 
tally are probably E c energy. 

It is possible to make the supposition that each of the 
examined energies executes some, expressly certain, 
function. For example, the main function of symmetric 
excitations can be activation of muscle proteins. At the 
same time, they can activate both membrane and en- 
zymatic proteins that are quite often actually observed 
in the activation of myosin [9-11]. 

Antisymmetric excitation energy is not enough to excite 
the muscle protein because it lies below the symmetric en- 
ergy. Activation of membrane proteins can be their main 



function. At the same time, these excitations are able to ac- 
tivate enzymatic proteins that are also actually observed 
often enough during activation of membranes [11-13]. 

And, lastly, asymmetrical excitations have only one 
function - to activate exceptionally enzymatic activity in 
those cases, when membrane and muscular activities are 
not needed. That is only for intracellular processes. 

Conformational response to the excitation of the 
alpha-helical region of the protein molecule 

For the analysis of conformational response of the 
alpha-helix on the considered excitations, it is necessary 
to appeal again to new equilibrium values of the step of 
the alpha-helix. From definition (3), it is possible to find 
Rna = Ro • (1 ~ fi\A a n\ 2 )> where designation is entered: 

B= ^fjyj : I . If we consistently apply the model of 

Rq • I W< 1 (Ro ) | 

dipole interaction between the peptide groups, then 
where, as mentioned above, Ad -0.29 D and d~3.7 D. 
Therefore, in this dipole model [14], ft ~ KT 1 . Taking into 
account the definitions of coefficients A aw given in (7), it 
is possible to get following: 

1. It is possible to obtain the following formula for 

symmetric excitations: = R 0 - (l— . That is, 

all three chains are reduced equally and evenly in 
the space. Then the length of every peptide chain 
can be appraised, so 

n—l 

This change is small and, at first glance, has no 
practical significance. But it will be so only in the 
classical model of the alpha-helix (Figure 2). If we 
consider, for example, that the peptide chains of 
myosin themselves form superhelices, then the 
effect of contraction increases. This is done by 
changing all characteristics of an alpha-helix: the 
step of the helix, its radius, and the effective 
number of peptide groups on the turn of the 
helix. Also, additional self-torsion takes place. The 
strengthening of the effect of contraction is deter- 
mined by the mutual torsion of long alpha-helical 
regions of light faction of myosin and their tor- 
sion on actin filaments. 

2. For antisymmetric excitations, it is possible to obtain 

*W = /£> = R 0 - (l- J^,R$ = R 0 . Respective 
lengths are as follows: 

T {a) _ T {a) _ T fiRo . T {a) _ t = n \t 

L Q —L x — Lq — , L 2 —Lq—KqN c . 
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In this type of excitation, one of the peptide chains 
does not change (here, it is a chain with the number 2), 
and two others are reduced up to the value 
Such asymmetry is enough for the alpha-helix to 
take a form of the segment of torus instead of 
cylinder (Figure 3). Application of the simple 
geometric considerations gives for the radius of 
curvature and angle cp: 



_2 d a _ 2 

Rk= ~o ' IT ' L ° = ~o ' d cc' N c', 
p K o p 



2 d a 



and for displacement A, it is possible to get such estimation: 



2V3 d a 2V3 d a 



(16) 



Taking into account the numerical values /5 ~ 10 1 , 
R 0 = 5.4 A, and d a = 4.56 A in (16) gives A- ^ (A). 
For the typical number of turns in many enzymes 
and membrane squirrel (N c > 10), displacement will 
have an order A > 2 A. This is consistent with the 
observed values [11]. 
3. For asymmetrical excitation, the following values are 

implemented: R%> = R% =R 0 -(l-^, R<£ = R 0 - 

(l— The corresponding lengths of peptide 
chains equal 



r(«) _ r(«) _ J j{h 



Lq- 



2fiR 0 



A Y 




X 



Figure 3 Explanation to estimation of displacement A of free 
(here upper) end of alpha-helix for antisymmetric excitations. 



The nature of the distribution of deformation along 
the peptide chain for this type of excitation is similar 
to that of the antisymmetric excitation. The only 
difference is that the chain, which in the previous 
case has not changed at all, now has shortening 
stronger than the other two. It is possible to 
estimate displacement for this case too: 



P ■ 1 



6N C 3N 2 r \R 0 



Here, A is the displacement for antisymmetric 
excitations, which is determined by Equation 16. 
Unlike displacement A, displacement A (w) 'directed' 
to the opposite side. Executing numerical estimates, 
it is possible to set that A (w) > A, if the number of 
turns in the alpha-helix N c < 14, but at N c > 14, we 
will have A (w) < A accordingly. 

Consequently, asymmetrical excitations demonstrate 
two very interesting features. First, it has the lowest 
energy and at diminishment of the number of turns 
N c , it falls down yet more. Second, a conformational 
response for this type of excitation is the biggest for 
N c < 14. This is typical for enzymatic proteins only. 

Conclusions 

The general methods [7,15-17] of description of the ex- 
cited states of the condensed environments were applied 
to the alpha-helix region of a protein molecule. The 
alpha-helix is considered as a nanotube, and excitations 
of the environment are described as quasiparticles. It is 
shown that three different types of excitation exist, and 
each of them is probably used by three different types of 
protein. The symmetrical type of excitation is used for 
muscle proteins, the antisymmetric type of excitation is 
used for membrane proteins, and the asymmetric type of 
excitation is used for enzymatic proteins. It is possible 
that some excitations of asymmetrical type exist, which 
are also used by enzymes. The estimations were done for 
displacements of the free end of the alpha-helix. The ob- 
tained displacements are in agreement with experimen- 
tal data. Therefore, the obtained results can be the basis 
of the interpretation of the functional properties of pro- 
teins characterizing their activity related to their con- 
formational changes [11]. 

Endnotes 

a Off-system unit of energy: 1 eV = 1.602 x 10" 19 J. 
b For example, in the myosin protein, the helical region 
has about 200 turns or up to 700 amino acids. 
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